Magnetic adatoms on graphene in the Kondo regime: an Anderson model treatment 
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We study theoretically the Kondo effect for a magnetic adatom on graphene using the Anderson 
model. Upon obtaining the Green's function of the impurity to higher order contributions in the 
hybridization, we calculated analytically the selfenergy in the presence of strong correlations. It is 
found that the Kondo resonance takes place in a narrow energy range of the impurity level around 
the Fermi energy which can be tuned by a gate voltage. We show that this range is linear in the 
Fermi energy and is significantly narrower than in the case for a normal metal. The origin of this 
behavior is traced back to the inherent properties of graphene, especially its linear dispersion. The 
singularity in the full Green's function is also analyzed with the help of a transparent geometrical 
method. The relations between the various selfenergies and the implications for the experimental 
observations are discussed . 

PACS numbers: 75.20.Hr, 73.22. Pr, 72.15.Qm, 81.05.ue, 71.20.Tx 



I. INTRODUCTION 

Recently, graphene, a two-dimensional monolayer of 
carbon atoms, has attracted much attention [ij. Its low 
energy band structure consists of two inequivalent Dirac 
cones located at the Brillouin zone corners 0. In these 
cones, the energy-dispersion relation is linear (relativis- 
tic), and the dynamics of the charge carriers is gov- 
erned by a massless Dirac-type equation. The physics 
of magnetic adatoms on graphene is under active re- 
search by means of the scanning tunneling spectroscopy 
first-principles calculation [(|, and many-body ap- 
proaches that deal with the Kondo effect [7Hl9{ . Theoret- 
ical studies on magnetic impurities coupled to 2D Dirac 
fermions may be initiated in the context of d-wave su- 
perconductors [2Cj-|22j]. New aspects are particular to a 
graphene host such as the possibility of doping via a gate 
voltage that moves the Fermi level fi away from the Dirac 
point [l6[ , and the different symmetries realized depend- 
ing on where the adatom is positioned on the graphene 
lattice [I?} ■ Several studies have dealt with the physics 
of magnetic impurities on graphene. A Hartree-Fock ap- 
proximation |8l4lC| based on the Anderson model was 
utilized to study the mean-field behavior which is only 
valid at T > Tk, where Tk is the Kondo temperature. 
The Anderson model with infinite Coulomb correlation 
(U) under a magnetic field was studied with the conclu- 
sion of a Fermi liquid behaviour [121 ]. On the other hand, 
the Kondo model (or sd model) is also used to address 
spin interactions between the impurity spin and the spins 
of the Dirac fermions in graphene. The anisotropic sin- 



gle channel Kondo model and a two-channel Kondo 
model [lH were both considered. However, for half spin, 
the former leads to a Fermi-liquid character; while the 
latter gives rise to a non-Fermi-liquid-like ground state. 
These two are essentially different [23]. Very recently, 
the Kondo model was employed to address the quantum 
criticality and an expression of the Kondo temperature 
was given (l9j . In Ref. [l7|, a detailed symmetry group 



cal model. The conclusion is that, different symmetries 
are realized depending on whether the adatom is above 
one carbon atom (case A) or on the center of the honey- 
comb (case B). In Ref. [la], we focused on the Anderson 
model for the case A and conducted a Slave-Boson (SB) 
calculation. The results showed that the full occupation 
survives in a broad gate voltage range until the Fermi 
level tuned by the applied gate voltage touches the edge 
of the impurity level, in which case the zero-field sus- 
ceptibility is finite. However, the mean-field SB method 
does not answer how the Kondo effect evolves with the 
temperature and other parameters. Therefore, to address 
these issues transparently, we calculate here the Green's 
function including higher orders in the hybridization by 
using the method developed in Ref. [24| and discuss how 
the graphene properties are reflected in the calculated 
physical quantities. 



II. THEORETICAL MODEL 

The tight-binding Hamiltonian of the electrons in 
graphene is given by 
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H.c. 
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where a\ a (a,i,cr) creates (annihilates) an electron with the 
spin a on the position of the sublattice A, a (bj,a) 
creates (annihilates) an electron with the spin a on the 
position Rj on the sublattice B, stands for summa- 
tion over nearest neighbors, and t is the nearest neighbor 
hopping energy. The Hamiltonian H g is rewritten in the 
momentum space as 



(q) 4 



(2) 



where 



-t J2 e^'" 5 * • The Hamiltonian © can be 



analysis was conducted to clarify the appropriate physi- easily diagonalized to give E±(q) — ±i|</>(q)|, which can 
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be linearized around the Dirac points in the Brillouin 
zone 

E ± (k) =±v F \k\, (3) 

where vp = 3ia/2 is the Fermi velocity of the electron 
(t ~ 2.3 eV [l|). Please, note k is now a momentum with 
respect to the Dirac points. 

When introducing the impurity, we consider the situ- 
ation where it is above one carbon atom of the graphene 
sheet. The symmetry for this case gives rise to a mixture 
of the two Dirac cones and the Hamiltonian reads 17 1 



E 



(4) 



where e k — hvpk, s stands for the valley index, fc c is 
a cut-off of the momentum so that D — hvpk c where 
D is the cut-off of the energy (the explicit value of D 
will be given below). The Hamiltonian of the impurity 
is taken as H im = e J2a ftU + T J2a n <y n 3i where f a 
is the annihilation operator; here the spin a is at the 
impurity site, and n a = f\fa- The total Hamiltonian is 



H 

is 



-He, where the hybridization Hamiltonian 



(5) 



with vq 



2tt 



fio is the unit cell area, and vq is the 



hybridization strength in energy dimension. 



A. The Green's function and the selfenergies 

1. The first order approach 

We define the retarded Green's function 
((A(t)\B(0))Y = -iB(t)({A(t),B(p)}) where 
A{B) are operators. ((A\B)) U + is the tu + transform 
of ((A(t)\B(0))) r , where w ± = w ± i0+. In the fol- 
lowing, we calculate the retarded Green's function and 
symbol the index uj + as u) so that we can obtain the 
retarded (advanced) Green's function when replacing u) 
by u in the derived Green's function. The standard 
equation of motion (2|| [26| of the Green's function 
reads u((A\B)) u = ({A,B}) + (([A,H]\B)) U . We 
then calculate the Green's function for the impurity 
t aa ,{u) = {{f a \fl,)) u by using 

u((fM)) u = ({UJl}) + ({[U,H]\ft>)U (6) 

where 

[U, H] = e / CT + Uf a n 9 + «o E / * , \4* T kv feV - ( 7 ) 

s' 

Thus, we find 

(w-co)WH = 5 aa , +U((f a n s \fl)) L0 

+ v Q J2 f dky/\k\((c sk *\fi,))^ (8) 
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FIG. 1: (color online) The real part of Eo(u; + ) 

Eq. jrrj. 

By employing the equation 

[c sfcCT , H] = e k c ska + v \/\k\U> 



1.0 



as given by 



we obtain 



((C.ka\fl)) u = 



Sk 



(9) 



(10) 



Therefore, the third term on the right hand side of Eq. 
reads Eo(<jj)t aa > (w), where the selfenergy is defined 



by E (w+) =jiEJdk\k\ 
calculated as [18| 



£fc) 1 and is readily 



-In- 



|£ 2 



wr|a;|0(.D- |w|) 



(11) 



E (w + ) = -7] 

where ri = §a "° . Let us discuss the structure of this 
selfenergy (see Fig. [IJ. The imaginary part is simply 
proportional to the absolute value of the energy with re- 
spect to the Dirac point, i.e. \uj\. The real part is shown 
in Fig. [TJ When to — > 0, the real part tends to zero as 
well which is a ty pical character in the marginal Fermi 
liquid theory [23]. When ui approaches the band edges, 
a singularity occurs. A similar singularity at the band 
edg e is met in the impurity problem for a normal metal 
[28( 1 . The difference is that the parity in energy space 
(with respect to the Fermi level that is assumed to be in 
the middle of the band) is even (odd) in a normal metal 
(graphene). Another observation is that the zero points 

occur at lj — and oj — ±^D. 

Having obtained the selfenergy we can write down the 
equation 

[u - e - EoHW = <W + UG{lo), (12) 

where G(u) = {{fan^f^,))^. We can further formulate 
the equation of G as 



(u — £i)G(cj) = {n&)5 aa > + vq 



dk 



S 

+ F 2 (cj) - F 3 (cj)}, (13) 
e + U, Fx{lj) = {(n<rc ska \fl,)) u , F 2 (u) = 



where e% 

((fUs ks fa\fl))u, and F 3 (lj) 
higher order Green's functions 



{{c\ ks hfMl,)) 



are 
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2. Accounting for higher order terms 

It is straightforward to obtain the following equations 
in the same way for the higher order Green's functions 



and 



(w - e k )F x = v ^\k\G + v J2[ dk'^/Wl x 

s' 

(<</k,^c sfeCT |/t,)) - {{c\, k ,-J s c ska \fl))) , (14) 



(us - e k )F 2 



V \/\k\G + {flc a ks)6a<r' 
s' 

{{c\ lk , s c sks Wl))) , (15) 



and 



(w - £ 2 + £k)F 3 = -v \/\k\G + (c\ ks f s )8 a<7 , 



+v a V / dk'y/W\ {{{c\ k - a c s , k , & U\fl,)) 



(16) 



where £2 = 2e + U is the energy for the double occupa- 
tion of the impurity. As usual, to obtain useful expres- 
sions we have to truncate the iterative equation by some 
app roximations . We follow the standard method in Refs. 
[24l l29l . 30] by assuming that the operator pair with the 
same spin indices can be pulled out of the Green's func- 
tion (GF) and regarded as an average. This assumption 
amounts to considering only spin conserved couplings 
within graphene or between graphene and the impurity. 
Physically, this amounts to saying that there is no mech- 
anism to flip the spin when it hops from the impurity to 
graphene or when it propagates in graphene. The equa- 
tion of motion method is able to capture the Kondo ef- 
fect only qualitatively but can unravel explicitly how the 
Kondo resonance develops in the local density of state 
[26l [3l| . We show one example of the decoupling scheme 
in the equation of F± , i.e. 

{{ftc s ' k > S Cska\.fl l ))u ~ (fsC s tWa)((c aka \fl,)) u . (17) 

The others can be done in the same way. After straight- 
forward manipulations, we can find the GF 



fn 



l + (U/£j)((n s )+A d 



u - e - S - (E/y«)(£o4r - E19) ' 
where lj = lj — e\ — £2, £2 = 2£o(w) — £o(£2 — 

(fa c sks) 



(18) 



dkJ\k 



( C skaf<?) 



to - e k 



LJ - £ 2 + £fc 



(19) 



E lff H = v 2 J2 [ dkdk'VW\VW\(^^ 



LJ - £ 2 + £fc 



(20) 



We employ the relation between the average and the 
GF 

(ab) = -±3 J du'f(u/){{b\a)) u ,+ , (21) 
where f(u>') is the Fermi function. Here we introduced 



,/± 



co' ± i0 + for the retarded (advanced) GF. It is 
easy to show {flc sks ) = (c\ ks f s ). And the GF of the 
graphene with an impurity has the form 



((Cska \ c s > k '<j'))ii 



u - £fe 



{lj - £ k )(u> - S k >) 



(22) 



In the calculation of the selfenergy Y,\ s we only keep the 
first term in Eq. (|22"]) . i.e. we use the GF of the pure 
graphene system without the impurity. In the following, 
we only focus on the limit U — > 00 in which case we infer 
the expression 



1 - (Ua) - A s 



lj - £ - (1 - A s )S - Sis ' 



(23) 



The selfenergy Ei ff can be calculated exactly (see Ap- 
pendix) 



Sia(w + ) = r, 



1. 



\D\ \^-D*\ (I w-n 



where ip is the Digamma function. 



3. Calculations of the spectral function A&{uj) 
An evaluation of the function Ag-(oj) yields 

M«>) = / dcj'f(iv')v 2 J2 J dk \ k \ 



(24) 



(lj - e k + t0+)(uj' - s k + i0+) 



(u - e k + iO+)(Lj> - e k - i0+) 
The denominator in the first term can be split into 



(25) 



1 



1 



1 



lj - e k + i0+ lj' - e k + i0+ 
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Because of the term to'— to in the denominator the integral 
is peaked at to' « w. However, the difference of the two 
terms in the square brackets will disappear in this case. 
Therefore, the whole contribution of this integral is less 
important than the second term in Eq. (|25[) and can be 
neglected. Splitting the product in the denominator of 
the second term in Eq. (|25[) . we find 

AS[U ] = ^ J W '- W -i0+ 

(26) 

In accordance with the consideration that the most con- 



tribution comes from 



to, we may simplify the ex- 



pression further by using 

S (w+) - £ (u/-) « -2wr)\u'\e(D - \to'\), (27) 



to obtain 

A s (lo+) 



— D 



dw'\u>W)t% s {u,'+) 
■J -to- i0+ 



(28) 



A similar relation was derived for a normal metal as a 
host material in Ref. [24| . 

At low temperatures, by using the same assumption as 
in Ref. (24|, we treat the GF as being a smoother part 
than the other term in the integrand. Thus, GF can be 
pulled out of the integral and it depends on to rather than 
on to' . Doing so yields 



A ff (o,+ )=t; ff ( W +)£ lff ( W +). 



(29) 



At high temperatures, the function A tends to be 
small. We assume a zero order Green's function, i.e. 
t ss (uj' + ) = (1 — (n a ))/(to' + — So), and substitute it into 
Eq. (|28p leading to the compact formula of the function 
A at high temperatures 



A _, +) = (l-(n ff ))[S lg (o;+)-S lg ( e +)] 

£o - W 



(30) 



III. DISCUSSIONS AND NUMERICAL 
CALCULATIONS 

A. The case of a vanishing A- function 

1. The Kondo resonance and the Kondo temperature 

From the calculation of the function A we concluded 
that it is related to the selfenergy Sis-. It is necessary 
to investigate the singularity of the A-function but we 
postpone it to the next subsection. From Eq. (|29|) . we 
know that the only possibility of having a singularity 
in A is via the GF t ss since the selfenergy Eig- is non- 
singular, except for T — > (and we do not consider the 
zero temperature). For a comparison, the corresponding 
selfenergy Sig- for a normal metal has a singularity when 
to — > fi at low temperatures. 
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FIG. 2: (color online) The digamma function is shown with 
different gate voltages and temperatures in (a) and (b), re- 
spectively. In (a), T = Kr 5 D, and r) = 0.02. In (b), /i = 0.LD 



and 77 = 0.02. A graphic method to solve Eq. (|32[) is shown 
in (c). The short-dash straight line is for the left hand side 
of Eq. (I32|) where eo = 0.1. The curves are for the right 
hand side of Eq. (132 p where the digamma function is calcu- 
lated numerically. The other parameters are rj = 0.02 and 
T = 10 _5 D. The energy unit is taken as D in all graphs. 



We note that the function A depends linearly on rj 
which is a small quantity for the graphene case. If the 
function A is not singular exactly, we can neglect it in the 
numerator and denominator since the terms with A are 
at least one order of magnitude less than the other terms 
(in rf). This approximation yields a zero order result 
and we infer, in the following, that it does not change 
qualitatively the buildup of the Kondo resonance. Thus, 
here we study the case of a vanishing A-function. The 
Green's function now reads 



to 



(31) 



<W(1 - (riff)) 



oj — £q — 1.5Sq — rj 



2 111 (2irT) 2 



— \ 
2 1 2-rriT I 



The resonances are determined by the positions where 
the real part vanishes, i.e. 



uj-e 



1.5RE + i] 



\D\ ]n \" 2 -D 2 \ 
2 (2ttT) 2 



(32) 



The Kondo resonance manifests itself in the solution 
around the Fermi level and is largely determined by the 
behavior of the digamma function. We thus study the 
behavior of digamma function at first. In Fig. [2a) and 
(b), the digamma function variation with to is shown for 
different gate voltages and temperatures. Compared with 
the results for a normal metal, we observe new features, 
i) The r\ is small in graphene while it is larger in normal 
metal because of the constant large DOS. ii) A factor 
appears in front of the digamma function which sup- 
presses the resonant effect when /1 tends to half-filling 
(fi — 0). This means the Kondo resonance is suppressed 
when the gate voltage is small and eventually disappears 
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at half-filling. This linear dependence on \n\ is just a 
consequence of the relativistic dispersion of graphene. A 
another pronounced effect of it is that the temperature 
dependence in the second term and the third term of the 
right hand side of the Eq. (1321) do not cancel with each 
other. The consequence is that there is no singularity 
when u> tends to the Fermi level. This is simply under- 
stood. For a very small but finite temperature, if uj — > fi 
then tp — > 0(1/2) rather than a logarithmic divergence in 
normal metal case (see Ref. [24J). This is reflected by the 
uniform topmost value of the peaks in Fig. [^a) and (b) , 
i.e. —'0(1/2) « 1.96351. The peaks are emerging at the 
Fermi level which can be tuned with the gate voltage (see 
Fig. [2ja)). The peaks only exist in a narrow range of cj in 
which \u — can be comparable to 2irT. It decreases dra- 
matically outside this interval. When \uj — fi\ ^> 2nT, the 
real part of the digamma function can be approximated 

by a logarithmic function, i.e. 3^-0 = In ^-kT^ 1 " 1 ■ 
In this case, the digamma function develops a singularity 
only when T — > 0. Please, note that this singularity in 
temperature is removed by the second part of Eq. Q32[) for 
a normal metal. However, the gate-voltage dependence 
of the digamma function lets this term survive (see Fig. 

EZtb)). 

Because of the smallness of r\ and the non-singular 
peak located at the Fermi level, it is relevant to inspect 
whether a solution around the Fermi level (i.e. the Kondo 
resonance) can be obtained. For a gate voltage \fi\ < D, 
3fJ£ does not give rise to a singularity. The second term 
in Eq. (|32p , with the logarithmic behavior increases with 
lowering the temperature rendering possible a Kondo res- 
onance. To show this clearly, we resort to a graphic 
method to solve Eq. (|3"2"j) in Fig. [He). Whenever an in- 
tersection of the peaks of the curves and the straight line 
occurs, the Kondo resonance emerges. We should em- 
phasize that the vanishing DOS of graphene at the Dirac 
points at half-filling suppresses the Kondo resonance. To 
release this suppression, a reasonable gate voltage is nec- 
essary. As a consequence of the non-singular peak at the 
Fermi level, we infer the important conclusion that the 
Kondo resonance can only be observed in a quite nar- 
row window of eo or the gate voltage with respect to the 
Fermi energy (at a fixed low temperature). 

Kondo temperature : If the temperature is very low and 
j3(u) — fj,) 3> 1, we may obtain a virtual level as 

e' = £o + 1.59fffio + r/M In ^— ±- . (33) 

V l £ o Ml / 

Since we seek a solution for the Kondo resonance near 
the Fermi level, we may from Eq. (|33[) obtain that 

|4 - « D D ^{2TTT)^ D ^e . (34) 

The estimation of the Kondo temperature in this case 
should be supplemented by studying the scaling process 
32]. In the scaling method, the bandwidth D is re- 
duced to an effective bandwidth D and the Hamiltonian 
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FIG. 3: (color online) LDOS, evaluated as DMj a , is shown 
for different eo in (a) through (f). fi = —0.1, T = 10 , and 
n = 0.02. £ = -0.45,-0.34,-0.32,-0.31,-0.28, and -0.1 
in (a) to (f) respectively. Energy is measured in D. n a is 
assumed to be 0.5. 

is renormalized with a new virtual level. This reduction 
can be done until D reaches the Fermi energy. In view of 
this picture, we take D — > D ~ \fj,\ and obtain the Kondo 
temperature Ik 

|4-/z|«He^f wT K . (35) 

The real part of Eo is ignored for the case where is 
not so large. There is no difficulty to include the Eo by 
setting its argument u> = ji. When tends to zero, the 
Kondo temperature tends to zero for eo < /i. This is a 
natural observation from Fig. Efc), and can be under- 
stood by the fact that no solution around the Fermi level 
can be obtained for a vanishing digamma function in this 
case. The intuitive picture is, there are no free electrons 
to screen the spin, since the local density of state is zero. 
With increa sing \/j,\, the Kondo temperature increases ex- 
ponentially [17] with a linear prefactor in This build 
up is qualitatively consistent with the prediction made 
on the basis of the sd model [13| . 

In Ref. , the authors considered 2D Dirac fermions 
with the linear DOS \oj\ t where r is a dimension in- 
dex (r = 1 for graphene). They studied the scaling 
and revealed that T +K (r -> 1) = \fi\e~ 1 - 1 '/^"^ and 
T_ K (r -4- 1) = \n\[W(l + r)/(l + r)} 1 ^, where ± for 
/i ^ 0. These results also show that the Kondo effect 
is sup pres sed when fi — > 0. By quoting the data from 
Ref. i-e- Tk ~ 15K at \fi\ = 0.2 eV, we esti- 

mate the impurity level to be eo ~ [i — 5ry|/i|. We note 
Tf\fi\ characterizes the hybridization strength at the Fermi 
level fj,. We first wish to compare this coupling with that 
characteristic for a normal metal or semiconductor quan- 
tum dot systems, i.e. A = itvq/(2D) [24|. The ratio is 
22T = niT6evp ' wri ere a = 1.42A and vf = 10 6 m/s have 
been used Q. If we take D « 7 eV 19] and \fx\ = 0.2 
eV, the ratio is about 0.0035. Therefore, we infer that 
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the coupling between graphene and the adatom impurity 
is weaker than in the case of a normal metal host. Sec- 
ondly, we estimate the value of 77. If vq = 1 eV (HI, we 
find 77 rj 0.0085. Thus, we estimate for the impurity level 
£0 = jtt— 8.5 meV (for \\x\ = 0.2 eV). The impurity level is 
hence quite close to the Fermi energy for an observation 
of the Kondo resonance. 




(a) (b) (c) 



B. LDOS of the impurity 

Let us define the local density of state (LDOS) for 
the spin a as Nf„{uj) — — ^t aa (uj + ), and show how the 
Kondo resonance develops and disappears when varying 
the impurity level in Fig. [31 We note that ignoring the 
function A in the numerator could induce a quantitative 
difference in the calculated occupation, even when A is 
small. However, in this case it will not induce a qualita- 
tive shift since it does not contribute significantly to the 
Kondo resonance. In the denominator, the term with the 
function A and S can be neglected since this is a second 
order term in rj, while 77 is much smaller than its coun- 
terpart in a normal metal. Under these considerations, 
we performed numerical calculations of the LDOS; the 
results are shown in Fig. |3l 

In Fig. GHa) through (f), the impurity level energy is 
increased gradually. In (a), the Kondo resonance is not 
yet developed. In (b) and (c), it starts emerging, and 
in (d) a pronounced Kondo peak has developed. In (e), 
when the virtual main peak crosses the Kondo resonance, 
the mutual interplay shows as a dip. In (f), the virtual 
main peak passes over the Fermi level. As a consequence, 
the Kondo resonance disappears. This build-up of the 
Kondo resonance shown in Fig. [3] is also valid if we view 
the increase of the impurity level energy as a lowering 
the Fermi level (while keeping fixed the impurity level 
energy). Therefore, our findings can be tested easily by 
tuning the gate voltage in the experiments. 

To gain more insights in the specific character of the 
Kondo resonance in graphene, we compare what we ob- 
serve here with the normal metal case where the Kondo 
temperature is T K ~ (DA) 1 / 2 exp[^S^i£)] m the U -> 00 
limit [24], H3- When /i is fixed, the Kondo temperature is 
decreasing when moving £q away from /i (below it). The 
Kondo resonance in LDOS disappears when Tk is below 
T, similar to the behaviour observed in Fig. [3] The dif- 
ference to the normal metal case is in the parameter win- 
dow for £0 where Kondo resonance exists. This window 
is for graphene rather narrow compared with the nor- 
mal metal, as the ratio ^ s quite small (about 0.0035, 
as given above). According to the parameters used in 
Fig. [3] and to keep T ^ Tk, the parameter window for 
£0 derived by using Eq. ([55)1 should be in a range of 
O77 1 I sa 0.02 below the Fermi level which is consistent 
with the observation from the figure. 



FIG. 4: (color online) A geometrical analysis for the Kondo 
resonance and the singularity of the Green function is shown. 
In (a) the function A vanishes. In (b) and (c) A is finite. The 
dashed vertical lines indicate the positions where the reso- 
nances can occur. R (I) on the axis means real (imaginary) 
axis. In (b) and (c), the rings are the singular rings with poles 
of the Green function. 



C. Non-zero A function study 

In this subsection, we wish to discuss the situation 
where the function A is solved by means of GF itself, 
i.e. Eq. (1291) . Thus, we find an iterative equation for t ss 
and the A^-function. We consider the spin degenerate 
case and infer 



t = — 



m — 



3 + t*£ Ei' 



(36) 



where 3 = z — Si, z = u) — eq — So, m = 1 — (n) and 
n = iv\- = 774. By substituting t*, we obtain the GF 



t = 



-A± VT 



2(SS >S S 1 )SJ' 
A = (|S| 2 -|S 1 | 2 ) + 2zm3(S S 1 ), 
T = (|S| 2 -|S 1 | 2 ) 2 + 4m(|S| 2 + |S 1 | 2 )3?(S S 1 ) 
- 4m 2 [3(S Si)] 2 - 8m|Si| 2 5R(SS;). (37) 

Please, note that T is a real function. The singularity 
occurs when 



SSq — SqSi — 0. 



(38) 



By using the conjugate relation 3* So = SqS| we obtain 
the constraint 



ISiI, 



(39) 



that leads to a singularity. To discuss this equation, let 
us first inspect the singularity analysis for a vanishing 
A. To do so we view a complex function as a vector in 
the complex plane (cf. Fig. [4]). The singularity in this 
case is given by |3| = or z = Si. That means, the 
vector z coincides with Si. However, when they are not 
completely coincident but merely have the same real part 
(which is shown by the two vectors reaching the vertical 
dashed line in Fig. 0Ja)) the vanishing real part gives 
rise to a prominent resonance. If it happens at the Fermi 
level, we observe a Kondo resonance. 
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For nonzero A, by noting the fact that 3 = z — Si, 
the relation in Eq. (|39p reflects a geometrical relation 
shown in Fig. EJb). A singular ring is generated with 
the radius set by the vector Si and we have 3 = z — Si. 
The basic physics in this singular ring is that the vector 
z and Si are varying with uj. At some particular values, 
the endpoint of the vector z lies on the ring defined by 
|Si|. A singularity is then encountered. Thus, the ring 
reflects the singularity of the GF in our investigation. 
Considering the point "O" on the ring, it means 5 = —Si 
or z — 0. Thus it reflects a resonance developed at the 
main peak of LDOS sitting at a renormalized virtual level 
e' = eo+5RSo. Let us consider the other special point, i.e. 
P, where 3 = Si or z — 2Si. Therefore, this point can 
be regarded as an extended point from the zero A case 
since the same singularity at the Fermi level is preserved 
(except for an enhancement factor 2). The picture is 
that introducing the A-function doubles the effect of the 
interaction and it is then easier to observe the Kondo 
resonance compared with the zero A case. 

So far, we discussed the meaning of the constraint con- 
dition in Eq. (|oT))) . Now let's discuss the meaning of Eq. 
([38]). If we define S = p e l \ we get 



3' = Ee~ 2ie 



Si, 



(40) 



where S', 3 and Si are viewed as vectors in the complex 
plane. Then the meaning of Eq. (PfUl) is that H' is ob- 
tained by rotating the vector 3 by an angle —29 while 
maintaining its magnitude. After this rotation, the new 
vector 3' (the "P" point on the ring in Fig. QJb)) should 
be equal to the vector Si. This situation is shown in Fig. 
S]Jb). When the end of the vector 3 is not lying on the 
ring for u — fi, say beyond it as shown in Fig. HJc), but 
the rotated vector 3' has the same real part of Si, we 
may observe a Kondo resonance. 
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Appendix A: The calculations of the selfenergy Ei 

We introduce 

"° de\e\f(e) 



$ ± (w) 



where ur 



I-d w - e 
to ± + . To calculate it, we define 

r-D 



de\e\[f(e)-l/2] 



-D 



The we make use of the relation 

= X ± M + (27 ? )- 1 So^ ± ). 



(Al) 



(A2) 



(A3) 



The calculation of \ are performed on an analytical con- 
tour which is shown in Fig. [5j 



H+3raTt Fermi 

I Energy 
H+niT | 



-D 



(I) 

I 

-D-2dT 



— ihf — ^ 

(i-iaTffl c K 11 ) 

: D-2niT 
H-37iiT? 



IV. SUMMARY 

In this work, we showed how the Kondo singlet de- 
velops for a magnetic adatom on graphene. We started 
from an Anderson model and kept the higher order con- 
tributions in the hybridization. Therefore, in this model, 
charge fluctuations are taken into account already. Ana- 
lytical Green's function is derived under a standard de- 
coupling scheme. The Kondo resonance can thus be stud- 
ied explicitly. It is found that the Kondo resonance takes 
place in a much narrow range for the impurity level with 
respect to the Fermi energy which can be tuned by the 
gate voltage. We show that this parameter window is 
proportional to the Fermi energy \fj,\ and is well narrower 
than that for the case of a normal metal host. Numeri- 
cal calculations for LDOS substantiate this analysis. The 
singularity in the full Green function is also analyzed in 
a transparent geometrical method, i.e. a singular ring 
in the complex plane. The relation between the various 
selfenergies is revealed. 



FIG. 5: (color online) The contour loop for the calculation 
of the selfenergy Eia is shown. The dashed path is deformed 
into the path indicated by solid lines with a circle C. 

Firstly, we need to clarify the meaning of |e| in an ana- 
lytical contour sense. From some further analysis we infer 
|3?(e)| when e is complex. The Fermi-Dirac distribution 
function has periodicity on the imaginary axis through 
2iriT (we set the Boltzmann factor ks = 1), where T is 
the temperature. Therefore, we find 



X ± {u + 2mT) 



D-2-kiT 



-D-2-iriT 



do/yi[/K)-i] 

lo-lo" + i0+ 



(A4) 



when we shift the argument by 2iriT. The integral path 
in Eq. (|A4[) can be deformed into —D - 2niT H> -D H> 
D D - 2niT + pole C. The pole is sitting at fi — iriT. 
Therefore, we get 

X + (oj + 2-KiT) = x + M + £ +(I) + (II), (A5) 
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where the second term is a circular integral around the Defining z — | + f^qy , Eq. 
pole along the path C, the third and the fourth terms are digamma function 
integrals along the path (I) and (II), respectively. The 
second term can be obtained by using the residue calculus 
as 



recovers a relation of the 



2niT 



i c ( w — m) + 7r *^ 1 

The third and the fourth terms can be calculated as 



(A6) 



(I) + (II) 



(A7) 



\D\ 1 /(oj- D + 2TriT + iO + )(u) + D + 2TriT + iO + ) 
~Y n V {uj - D + iO+)(uj + D + i0+) 

Thus, we obtain the relation 

$+ (uj + 2niT) \D\ 

(u>- D + 2niT + iO+)(uj + D + 2niT + i0+) 
(2mTf 

$+(w) \D\ , / (uj - D + iO+)(uj + D + i0+) 
In 1 



x In 



2mT 
(uj — fi) + iriT 



(2niTy 



(A8) 



if>(z + l)=il>(z) + -, 



(A9) 



leading to 



ir>l 



In 



|w 2 - D 2 



(2irTf 



1 , u_ - M 

2 2?riT ' 



(A10) 



From the definition of Ei^(w ± ), we obtain 

s lff ( w ± ) = 7 7 $ ± H. 



(All) 



Substituting the corresponding expressions, we find Eq. 

(HI. 
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